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Abstract: The basic work of Zaslavskii et al showed that the classical non- relativistic 
electromagnetically kicked oscillator can be cast into the form of an iterative map on 
the phase space; the resulting evolution contains a stochastic flow to unbounded energy. 
Subsequent studies have formulated the problem in terms of a relativistic charged particle 
in interaction with the electromagnetic field. We review the structure of the covariant 
Lorentz force used to study this problem. We show that the Lorentz force equation can be 
derived as well from the manifestly covariant mechanics of Stueckelberg in the presence of 
a standard Maxwell field, establishing a connection between these equations and mass shell 
constraints. We argue that these relativistic generalizations of the problem are intrinsically 
inaccurate due to an inconsistency in the structure of the relativistic Lorentz force, and 
show that a reformulation of the relativistic problem, permitting variations (classically) in 
both the particle mass and the effective "mass" of the interacting electromagnetic field, 
provides a consistent system of classical equations for describing such processes. 
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1. INTRODUCTION 

Zaslavskii et al [1] have studied the behavior of particles in the field of a wave packet 
of an electric field in the presence of a static magnetic field. For a broad wave packet 
with sufficiently uniform spectrum, the problem can be stated in terms of an electrically 
kicked harmonic oscillator. They find that for rational ratios between the frequency of the 
kicking field and the Larmor frequency associated with the magnetic field, the phase space 
of the system is covered by a mesh of finite thickness; inside the filaments of the mesh, the 
dynamics of the particle is stochastic and outside (in the cells of stability), the dynamics 
is regular. This structure is called a stochastic web. It was found that this pattern covers 
the entire phase plane, permitting the particle to diffuse arbitrarily far into the region of 
high energies (a process analogous to Arnol'd diffusion [2]). 

Since the stochastic web leads to unbounded energies, several authors have considered 
the corresponding relativistic problem. Longcope and Sudan [3] studied this system (in 
effectively 1 + 1/2 dimensions) and found that for initial conditions close to the origin of 
the phase space there is a stochastic web, which is bounded in energy, of a form quite 
similar, in the neighborhood of the origin, to the non-relativistic case treated by Zaslavskii 
et al [1]. Karimabadi and Angelopoulos [4] studied the case of an obliquely propagating 
wave, and showed that under certain conditions, particles can be accelerated to unlimited 
energy through an Arnol'd diffusion in two dimensions. 

The equations used by Longcope and Sudan[3] and Karimabadi and Angelopoulos [4] 
are derived from the well-known covariant Lorentz force 

'"=•» isH^f- (L1) 

where ds is usually taken to be the "proper time" of the particle. Multiplying both side 
by dx^/ds and summing over \i (we use the Einstein summation convention that adjacent 
indices are summed unless otherwise indicated, and the metric is taken to be (— , +, +, +) 
for the indices (0, 1, 2, 3), distinguishing upper and lower indices), one obtains 

dXfj, d 2 x fI _ 1 d , dxj, dx^ . _ 

ds ds 2 ~ 2 ds [ ds ds ] ~ ' 1 ' ' 

taking the usual value for the constant (in s), we have that 

1T1T = < L3 > 

This result provides a consistent identification of the parameter s on the particle trajectory 
(world-line) as the "proper time" : 

ds 2 = — Ttdx^dx^ 
a 1 j,,2 



= dt z - -^dx 2 (1.4) 
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so that 

ds 

dt — — ^=^= = 'yds, (1-5) 

the Lorentz transformation of the time interval for a particle at rest to the interval ob- 
served in a moving frame. This formula has been used almost universally in calculations of 
the dynamics of relativistic charged particles [6,7]. The Lorentz trtansformation, however, 
applies only to inertial frames. Phenomena occurring in two inertial frames in relative mo- 
tion are, according to the theory of special relativity, related by a Lorentz transformation. 
An accelerating frame, as pointed out by Landau and Lifshitz[6], induces a more compli- 
cated form of metric than the flat space (—, + , + ,+). Mashoon [11] has emphasized that 
the use of a sequence of instantaneous inertial frames, as has also often been done, is not 
equivalent to an accelerating frame. He cites the example for which a charged particle at 
rest in an inertial frame does not radiate, while a similar particle at rest in an accelerating 
frame does. As another example, consider again the first of (1.4). If we transform to the 
inertial frame of a particle with constant acceleration along the x direction, 

/ 1 2 

x = x + -at , 

then (1.4) becomes (as in the discussion of rotating frames in [8]) 

12 1 

ds 2 = (1 - -^a 2 t 2 )dt 2 - -^atdx'dt - ~^{dx' 2 + dy 2 ). 

In the frame in which dx' = dy = dz = 0, dt is the interval of proper time, and it is not 
equal to ds. For short times, or small acceleration, the effect is small. We shall discuss 
this problem further in Section 3. 

Continuing for now in the standard framework, Eq. (1.3) effectively eliminates one of 
the equations of (1.1). We may write 

<i 2 x , dt x d ,dtdx., 

ds 2 ~ ^ds'Jt^TsW 

e dt 1 „, 
= E+ -v x H). 

m ds c 

Cancelling 7 = dt/ds from both sides, one obtains 

i 7 v) = -(E+-vxH), (1.6) 

dt m c 

the starting point for the analysis of Longcope and Sudan [3] and Karamabadi and An- 
gelopoulis [4]. A discrete map can be constructed from (1.6) just as was done for the 
nonrelativistic equations of Zaslavskii et el [1]. As we have remarked above, the stochastic 
web is found at low energies; it deteriorates at high energies due to the 7 factor. 
The time component of (1.1) is 

d 2 t _ e E dt ^ 

ds 2 mc ds 
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or 

(1.8) 

at mc z 

Landau and Lifshitz [6] comment that this is a reasonable result, since the "energy" of 
the particle is 7mc 2 , and eE • v is the work done on the particle by the field. It is 
important for what we have to say in the following that Eq. (1.7) is not interpretable in 
terms of the geometry of Lorentz transformations. The second derivative corresponds to 
an acceleration of the observed time variable relative to the "proper time" ; the Lorentz 
transformaton affects only the first derivative, as in (1.4). We understand this equation 
as an indication that the observed time emerges as a dynamical variable. Mendonga and 
Oliveira e Silva [8] have studied the relativistic kicked oscillator by introducing a "super 
Hamiltonian" , resulting in a symplectic mechanics of Hamiltonian form, which recognizes 
that the variables t and E are dynamical variables of the same type as x and p. This 
manifestly covariant formulation is equivalent to that of Stueckelberg [9] and Horwitz and 
Piron [10], which we shall discuss in the next section. 

We have computed solutions to the Lorentz force equation for the case of the kicked 
oscillator (see fig. 1), using methods slightly different from that of Longcope and Sudan 
[3] and Karimabadi and Angelopoulis [4]. At low velocities, the stochastic web found by 
Zaslovskii et al [1] occurs; the system diffuses in the stochastic region to unbounded energy, 
as found by Karimabadi and Angelopoulis [4]. The velocity of the particle is light speed 
limited by the dynamical equations, in particular, by the suppression of the action of the 
electric field at velocities approaching the velocity of light [5] . 

The rapid acceleration of the charged particles of the kicked oscillator further suggest 
that radiation can be an important correction to the motion. The counterexample of 
Mashoon [11] was based on the phenomenon of radiation. It has been shown by Abraham 
[12], Dirac [13], Rohrlich [6] and Sokolov and Ternov [14] that the relativistic Lorentz force 
equations in the presence of radiation reaction is given by the Lorentz- Abraham equation 



mx » = -F» v x v + -—mi—x^ - ^x»x u x v ), (1.9) 
c 3 c as c 2 



where ro = e 2 /mc 2 , the classical electron radius, and the dots refer here, as in (1.1), to 
derivatives with respect to s. Note that from the identity (1.3), it follows (by differentiation 
with respect to s) that 

x^ = 0, + = °» ( L1 °) 

and hence (1.9) can be written as 

mx » = -F^x u + -^?tmW + 4rX M iv). (1.11) 
c 3 c as c z 

The last factor on the right is a projection orthogonal to x^ (if = — c 2 ), and therefore 
(1.11) is consistent with conservation of x^x^. Sokolov and Ternov [4] state that this 
conservation law follows automatically from (1.9), but it is apparently only consistent. 
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Radiation reaction therefore also implies that the connection between proper time and the 
Lorentz invariant interval may be subject to question. 

We have calculated the motion of the kicked oscillator using the form (1.9) of the 
Lorentz force, corrected for radiaton reaction, undoubtedly a good approximation under 
certain conditions, and will report on this in another paper in this volume [15]. 



2. THE STUECKELBERG FORMULATION 

As we have remarked above, Mendonga and Oliveira e Silva [8] have used a "super 
Hamiltonian" formulation to control the covariance of the electromagnetically kicked oscil- 
lator. Their formulation of the problem is equivalent to the theory of Stueckelberg [9] and 
Horwitz and Piron [10]; we shall therefore use the notation of the latter formulation. We 
first explain the physical basis of this theory, and then derive sytematically the covariant 
Lorentz force from a model Hamiltonian. 

The original thought experiment of Einstein [16] discussed the generation of a sequence 
of signals in a frame F, according to a clock imbedded in that frame, and their detection 
by apparatus in a second frame F' in uniform motion with respect to the first. The time of 
arrival of the signals in F' must be recorded with a clock of the same construction, or there 
would be no basis for comparison of the intervals between signals sent and those received. 
It is essential to understand that the clocks in both F and F' run at the same rate. The 
relation of the interval Ar between pulses emitted in F and the interval between signals 
At' received in F', according to the (equivalent) clock in F' is, from the special theory of 
relativity, given by 



This time interval, measured on a "standard" time scale established by these equivalent 
clocks, is identified to the interval At', the time interval between signals in the first frame, 
observed in the second, and called simply the time by Einstein. One sees that this Ein- 
stein time is subject to distortion due to motion. In general relativity, it is subjuct to 
distortion due to the gravitational field as well, and in this case the distortion is called the 
gravitational red-shift. We see that there are essentially two types of time, one correspond- 
ing to the time intervals at which signals are emitted, and the second, according to the 
time intervals for which they are detected. The first type, associated with signals that are 
pre-programmed, is not a dynamical variable, but a given sequence (as for the Newtonian 
time) , and the second, associated with the time at which signals are observed (the Einstein 
time), is to be understood as a dynamical variable both in classical and quantum theories 
[17]. 

Stueckelberg [9] noted that for a free particle, the signals emitted at regular inter- 
vals would be recorded at regular intervals in a laboratory, since the free particle would 
be in motion with respect to the laboratory with the same relation as between F and 
F'; the motion would then be recorded as a straight line (within the light cone) on a 
spacetime diagram. In the presence of forces, however, this line could curve. A sufficient 
deviation from the straight line could make it begin to go backward in time, and then 
the coordinate t wuld no longer be adequate to parametrize the motion. He therefore 
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introduced an invariant parameter r along the curve, so that there would be a one-to-one 
corrrespondence between this parameter and the spacetime coordinates. He proposed a 
Hamiltonian for a free particle of the form (the parameter M provides a dimensional scale, 
for example, in (2.5); it may also be considered as the Galilean target mass for the variable 
(l/c)v/£ 2 -c 2 p 2 ) 

K = V ^ (2.2) 
2M v ' 



for which the Hamilton equations (generalized) give 

dx» _ OK 

~dr ~ dp^ ~ ~M 



(2.3) 



It is clear that such a theory is intrinsically "off-shell" ; the variables p and p° = E/c are 
independent, as are the observables x and t, so that the phase space is eight-dimensional. 
Dividing the equation for the space indices by the equation for the time index, one obtains 

v =«= c l- (2A > 

precisely the Einstein formula for velocity. Furthermore, for the time component, 

* E (2-5) 



dr Mc 2 ' 



in case the particle is "on- shell", so that Mc 2 = ^E 2 - c 2 p 2 , (2.5) reads 

dt 1 



coinciding with (2.1). 

Stueckelberg [9] then considered adding a potential term V(x), to treat one-body 
mechanics, and the gauge substitution p^ — > p fJ, —eA fJ, (x) for the treatment of problems with 
electromagnetic interaction. He proposed a quantum theory, for which the Hamiltonian 
generates a Schrodinger type evolution 

i%^ T {x) = K^ T {x). (2.6) 

Horwitz and Piron [10] generalized the Stueckelberg theory for application to many 
body problems. They assumed that the standard clocks constitute a universal time, as for 
the Robertson- Walker time (the Hubble time) of general relativity [18], so that separate 
subsystems are correlated in this time. In this framework, it became possible to solve, for 
example, the two body problem in both classical [10] and quantum theory [19]. 

The equations (1.1) are not generally derived rigorously from a well-defined Lagrangian 
or Hamiltonian. They result from a relativistic generalization of the nonrelativistic Lorentz 
force (which is derivable from a nonrelativistic Hamiltonian). In the following, we shall 
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derive these equations rigorously from the Stueckelberg theory, to emphasize more strongly 
the nature of the problem we have discussed above, and to clarify some important points. 

The Hamiltonianian form for a particle with electromagnetic interaction proposed by 
Stueckelberg [9] is 

K = 2M • (2 - 7) 

The equation of motion for x^ is (we use the upper dot from now on to denote differentiation 
by r, the universal invariant time) 

= 8K = (p^-f^aO) 

dp^ M K 1 

and we see that then 

dx'dXy _ 2( ds^2_ (^-f^(x))(p M -f^(x)) 

dr dr ~ ° [ dr> ~ M ' [ ^ } 

a quantity proportional to K, and therefore strictly conserved. In fact, this quantity is the 
gauge invariant mass-squared: 

(pM _ lA"(x))(pp - - c A^x)) = -m 2 c\ (2.10) 

where we define m as the dynamical mass, a constant of the motion. It then follows that 

2 ,ds^2 o/dt,2 /dx N 2 m 2 c 2 

<%> = c %) = ~Kp~ (2 ' n) 



and, extracting a factor of (dt/dr), 



(#) 2 = ^. (2.12) 



Up to a constant factor, the Stueckelberg theory therefore maintains the identity (1.3). 

We now derive the Lorentz force from the Hamilton equation (this derivation has also 
been carried out by C. Piron [20]). The Hamilton equations for energy momentum are 

dp^_ = _dK_ = {p v '-\A V ) ,edA 
dr dx^ M l c5x/ 

edx»dA v [ } 



c dr dx^ 



Since p M = M^- + -A^, the left hand side is (A M is evaluated on the particle world line 

s"(r)) 

dp" ^d 2 x» edA^dx u , n ^ AS 

lF = M ^ + ^!7' (2 ' 14) 
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and hence 



or 



where (d^ = d/dx^) 



M d 2 x^ _ e,dA v _ dA^.dx" 
dr 2 c dx^ dx v dr 

,A 2 x» e^ a dx v , rt . 



F i*v = d»A u -d u A». (2.16) 



The form of (2.15) is identical to that of (1.4), but the temporal derivative is not with 
respect to the variable s, the Minkowski distance along the particle trajectory, but with 
respect to the universal evolution parameter r. 

One might argue that these should be equal, or at least proportional by a constant, 
since the proper time is equal to the time which may be read on a clock on the particle in 
its rest frame. For an accelerating particle, however, one cannot transform by a Lorentz 
transformation, other than instantaneously, to the particle rest frame. It appears, there- 
fore, that the formula (2.15) could have a more reliable interpretation. The parameter of 
evolution r does not require a Lorentz transformation to achieve its meaning. 

Since m 2 is absolutely conserved by the Hamiltonian model (2.7), however, we have 
the constant relation 

m 

ds = —dr, (2.17) 

assuming the positive root (as we shall also do for the root of (2.12); we do not wish to 
discuss the antiparticle solutions here). Eq. (2.15) can therefore be written exactly as 
(1.1). 

We see that the Stueckelberg formulation in terms of an absolute time does not avoid 
the serious problem of consistency that we have pointed out before. It is clear that the 
difficulty is associated with the fact that the Stueckelberg Hamiltonian, as we have written 
it, preserves the mass-shell, and we therefore understand the identity (1.3) as a mass-shell 
relation. 

Returning to the Stueckelberg-Schrodinger equation (2.6), we see that the gauge in- 
variant replacement p^ — > p^ — |A M (x) is not adequate. The additional derivative on the 
left hand side of (2.6) must also be replaced by a gauge covariant term. The possibility 
of t dependence in the gauge transformation implies that the gauge fields themselves may 
depend on r. The gauge covariant equation should then be [21] 

*§pMx) = {^(Pm - - - 7«5}Vv0r), (2.18) 

, where the fields a a , ol = (0, 1, 2, 3, 5) with d$ = d/dr, change under the gauge transfor- 
mation ip — > expi^-Aif; according to a a — > a a — d a A. It follows from this equation, in a 
way analogous to the Schrodinger non-relativistic theory, that there is a current 

% = - oXrWP" - i-^Wr ~ Vv(^ + i-a^Wr}, (2-19) 
ZM c c 
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which, with 
satisfies 

d T p T +d^ = d a r=0. (2.20) 
We see that for p T — > pointwise (J p T (x)d 4 x = 1 for any r), 



J^(x) = / ft{x)dT (2.21) 

•/ — oo 

satisfies 

d^{x) = 0, (2.22) 

and can be a source for the standard Maxwell fields. Since the field equations are linear, 
with source j a , one identifies the integral J dra fJ '(x, r) (or, alternatively, the 0-mode) with 
the Maxwell potentials [21]. It then follows that the so-called pre-Maxwell fields a a have 
dimension L -2 , and that the charge eo has dimension L. The Lagrangian density for the 
fields, quadratic in the field strengths (a,/3 = 0, 1, 2, 3, 5) 

jaf3 = Qa a f3 _ Qp a a ^ 

which has dimension L -3 , must carry a dimensional parameter, say A, and from the field 
equations Xd a fa = e oi^, one sees that the Maxwell charge is e = eo/A [21]. 

We understand the operator on the right hand side of (2.18) as the quantum form of 
a classical evolution function 

K = ^(p„ - % M )(p" - - % 5 . (2.23) 

It follows from the Hamilton equations that 

£E_ = p ; a (2.24) 
dr M v ' 

and 

dp^ eo dx u da v eo da$ 
dr c dr dx^ c dx^ 

Hence, 

<Px» _e Q dx v e 0( da 5 da*. 

M H^~ T~fr u + 7^ " 77 j - [Llb) 

If we define x 5 = r, the last term can be written as d^a^ — d^a^ = f§ , so that 



Note that in this equation, the last term appears in the place of the radiation correction 
terms of (1.9). It plays the role of a generalized electric field. Furthermore, we see that 
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the relation (1.3), consistent with the standard Maxwell theory, no longer holds as an 
identity; the Stueckelberg form of this result (2.11) in the presence of standard Maxwell 
fields, where m 2 is conserved, is also not generally valid. We now have 



and does not vanish. The right hand side corresponds to mass transfer from the field to 
the particle. 

As for the method of Longcope and Sudan [3], we may transform the derivatives with 
respect to r to derivatives with respect to t in the equation (2.26) as follows. Defining 
C = dt/dr, it follows from (2.26) that there is an additional term in the analogous form of 
the rate of change of ( (we use lower case to denote the pre-Maxwell field strengths) , 

The space components of (2.26) can be written as 

d 2 x j eo r ; 1 



*» - C m[ £ + c (vxh) 



(2.29) 



To illustrate some of the properties of this system of equations, we treat a simple 
example in Appendix A. The effective additonal forces include not only the term associated 
with the work done by the field, but additonal terms associated specifically with the r- 
dependence of the fields, and the fifth (scalar) field 05. Given the fields f a p, Eqs.(2.28) 
and (2.29) form a nonlinear coupled system of equations for the particle motion. 

For a gauge (generalized Lorentz) in which d a a a = 0, the field equations [21] 

d*f a = ef 



become 

-d a d a a? = ef, 

where, classically, j 13 = x fJ, 5 4 (x — x(t)), p = S 4 (x — x(t)), and x = x^(t) is the world line. 
The analysis of these equations is in progress. 

It has recently been shown that, with the help of the Green's functions for the wave 
equations of the fields in x M ,r, that the self-reaction derived from the contributions on 
the right hand side of (2.26) are precisely of the form of the radiation reaction terms in 
the Abraham-Lorentz equations (1.9) in the limit that the theory is constrained to mass 
shell, i.e., that (1.3) is enforced [22]. The off-shell corrections provided by (2.26) make 
the system of equations consistent, and should therefore provide a basis for computing 
problems involving the interaction of radiation with relativistic particles in a consistent 
way. 
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CONCLUSIONS 



We have shown that the standard relativistic Lorentz force equations are not consistent 
since they imply the mass-shell constraint x^x^ = — c 2 , a relation that can be valid only 
for a charged particle moving at constant velocity. The corrections are generally small for 
short times or small accelerations, and therefore calculations made with this Lorentz force 
are, in many applications, quite acceptable. However, for very large accelerations (at large 
compared to c), they could become inaccurate. 

A consistent theory may be constructed from a fully gauge covariant form of the 
Stueckelberg [9] [10] manifestly covariant dynamics, a theory which introduces a fifth gauge 
field [21]. The Lorentz invariant force equation derived from this theory contains an addi- 
tional term which enters in a way analogous to the radiation reaction term in the Abraham- 
Lorentz-Dirac equation (the self reaction force derived from this generalized equation in 
the mass-shell limit coincides with the radiation reaction term obtained by quite differ- 
ent methods for the Abraham-Lorentz-Dirac equation; it contains contributions from both 
terms on the right hand side [22]). 

It appears that the consistency of the classical equations governing the interaction of 
charged particles with electromagnetic radiation requires that both the particles and the 
fields must be permitted to move "off-shell" , as in the vertices of quantum field theory. 
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Appendix A 

The purpose of the following example is to show that in some cases the fifth field / M 5 
can cause to an effect which is very similar to the radiation effect that is calculated by 
Lorentz-Dirac equation. The fact that the mass is not conserved (the off-mass-shell case) 
is equivalent, in the case of radiation, to loss of energy through the radiation process. The 
particular example that we treat here is that of a charged particle in the presence of an 
uniform magnetic field in z direction (B = (0, 0, Bo)). 

As for the radiation reaction term of the Lorentz-Dirac equation, we choose the fifth 
field term to be * 

r 5 = (Cii,C 2 ±,C 2 y,0), (A.1) 

where the dot indicates derivative with respect to r. The Lorentz force (2.26) can be 
written as a set of differential equations, 

ll d 2 te^dt . t . 



* It appears that for the usual form of the radiation reaction, in an example with the 
field magnitudes that we shall choose, the ^f— term seems to be negligible, and the x^'x^ 
may be approximated by a constant number; one is left with the x^ term. We choose the 
fifth field term to have a similar structure. This choice is appropriate due to the close 
relation of these the radiation reaction terms of the usual theory [22]. 
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The solution of Eq. (A2) is 



, , d 2 x eB dy e dx , . „. 

M ^ = -r^ + i c ^Tr < A3 > 

llT d 2 y eB dx e ^ dy .. 



i = i e (XlT , (A5) 



where a\ = Using the complex coordinate u = x + iy, Eqs. (A 3) and (A. 4) can be 
written as 

du 

— = — lilU + a 2 U, {A.O) 

dr 

where a 2 = and O = (the Larmor frequency). The solution is 

u = u exp a2T e~ mT . (A.7) 

Using u(r) one finds that, 



x = e a2T (xo cos(fir) + yo sin(Or)) 
y = e a 2 T (_^ sin(Or) + yo cos(Or)). 



(A8) 



As expected, the radiation effect is determined by the constants ol\ and a 2 . 

It is possible to calculate the actual velocities by dividing Eqs. (A.8) by t; this results 

in 

dx _ aT ( ( dx\ , f dy s 



*=^lUJ cos(SiT)+ UJ„ sin(Sir: . 

l--"(-(l) ^ + (l) ^ 

where a = ol\ — a 2 . Notice that when ol\ = a 2 , there is apparent radiation (decrease of 
amplitude) as a function of r but not as a function of t ; in terms of t (which is redshifted) 
the particle appears to be circling forever on the same circle. This remarkable illustration 
is somewhat analogous to the phenomenon in which there is an infinite time required for 
a particle to arrive at the Schwarzschild radius in the Schwarzschild coordinate t, but a 
finite interval in the proper time of the particle. 

The magnitude of the (t-) velocity of the particle is 

v = v e- aT . (A10) 

When a = — , where tq = —^2 (70 is the radiation constant of the Lorentz-Dirac equation), 
Eq. (A 10) is exactly the solution which was obtained using the Lorentz-Dirac equation 
[14] [15]. This result is consistent with the approximations we have made in constructing 
the example. 
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Figure Caption 

fig. 1. A typical relativstic stochastic web. 
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